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Abstract 

This article is the complement to [7], which proves that flows (as introduced by [5]) can be found efficiently 
for patterns in the one-way measurement model which have non-empty input and output subsystems of the 
same size. This article presents a complete algorithm for finding flows, and a proof of its' correctness, without 
assuming any knowledge of graph-theoretic algorithms on the part of the reader. This article is a revised 
version of [4], where the results of [7] also flrst appeared. 

1 Introduction 

In the one-way measurement model [1, 2, 3], algorithms are essentially described by a sequence of single-qubit 
measurements (where the choice of measurement may depend on earlier measurement results in a straightforward 
way) performed on a many-qubit entangled state. This many-qubit state may be described in terms of the state 
of an input system I , together with a graph G of entangling operations involving / and a collection of auxiliary 
qubits prepared in the |-|-) state: each edge of G represents a single controlled-Z operation between two qubits. 
After the sequence of measurements, any qubits left unmeasured still support a quantum state, and are interpreted 
as an output system O . A triple (G, /, O) belonging to a given pattern is called the geometry of the pattern. 

In [4], it was shown that the flow property defined by Danos and Kashefi [5] can be efficiently tested for a 
geometry (G,/, O) when |/| = \0\ . The property is the existence of a causal flow^, which describes a partial 
order ^ describing an order (independent of measurement angles) in which the qubits of the geometry may be 
measured to perform a unitary embedding, once suitable corrections are applied to the output qubits. Causal 
flows may allow quantum algorithms to be devised in the one-way measurement model without using the circuit 
model: [6] proposes one way in which this might be done. 

This article presents a complete algorithm for finding causal flows for a geometry (G, /, O) with |/| = \0\ in time 
0{km) , where k = \I\ = \0\ and to — \E{G)\ , suitable for an audience with no experience in graph-theoretic 
algorithms. This is a revised version of [4], re- written with the aim of focusing on the algorithm for finding flows 
for the sake of reference. For the graph-theoretical characterization of flows, this article refers to [7], which is an 
improved presentation of the graph-theoretic results presented originally in [4] . 

Although no knowledge of graph-theoretic algorithms is assumed, a basic understanding of graph theory and the 
one-way measurement model is essential. For basic definitions in graph theory, readers may refer to Diestel's 
excellent text [9]; I will use the conventions of [7, 8] for describing patterns in the one-way model. 



2 Preliminaries 

In this section, we will fix our conventions and review the results and terminology of [7]. 

*This work was supported in part by ARDA, ORDCF, MITACS, and CIAR. 

-"^ These arc simply called "flows" in [5]: I use the term "causal flow" in this article to maintain consistency with [7]. 
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2.1 Basic notation and conventions 

For a graph G , we write V{G) for the set of vertices and E{G) for the set of edges of G . Similarly, for a directed 
graph (or digraph) D , we write V{D) for the set of vertices and A{D) for the set of directed edges (or arcs) of 
D . li X and y are adjacent, we let xy denote the edge between them in a graph, and x ^ y denote an arc from 
X to y in a digraph. We will use the convention that digraphs may contain loops on a single vertex and multiple 
edges between two vertices, but that graphs cannot have either. 

When a graph G is clear from context, we will write a; ~ y when x and y axe adjacent in G , and write S"^ to 
represent the complement of a set of vertices S C V{G) . 

If C is a collection of directed paths (or dipaths), we will say that a; — > y is an arc ofC, and that the edge xy is 
covered by C , when a; — > y is an arc in a path P & C . 

In this paper, N denotes the non-negative integers. For any n G N , [n] denotes the set { j e N \ j < n} . 

2.2 Results for Causal Flows 

2.2.1 Definition and motivation 

Definition 1. A geometry {G,I,0) is a graph G together with subsets I,0 C V{G) . We call / the input 
vertices and O the output vertices of the geometry. A causal flow on (G, /, O) is an ordered pair (/, , with a 
function / : O"^ — > /"^ and a partial order ^ on V{G) , such that 

(Fz) a;~/(x); (Fii) x^f{x); (Fm) y ^ f{x) =^ x =4 y , (1) 

hold for all vertices x G and y G V{G) . We will refer to / as the successor function of the causal flow, and 
=^ as the causal order of the causal flow. 

A geometry (G, /, O) represents the information of a one-way measurement pattern which is independent of the 
order of operations and measurement angles. G is the entanglement graph of the pattern, / is the set of qubits 
which are not prepared in a fixed state initially (their joint initial state in the algorithm may be arbitrary), and 
O represents the set of qubits which are not measured in the pattern (which thus support a final quantum state). 

The conditions (Fi) - (fiii) arc motivated by how byproduct operators and signal dependencies arc indiiccd by 
commuting correction operations to the end of a pattern which performs a unitary embedding. The significance 
of a causal flow on a geometry (G, I, O) is that any pattern defined on that geometry can be transformed into one 
which has the same measurement angles and which performs a unitary embedding Tij — > 7io ■ In particular, 
this means that unitary embeddings can be devised in the measurement model by ignoring signal dependencies 
and treating each measurement operator as though it post-selects for some one of the states in the basis of the 
measurement. See Section 2.2 of [7] for details. 

2.2.2 Grapii-theoretic characterization 

The result of [7] was obtained by characterizing causal flows in terms of collections of vertex-disjoint paths. 

Definition 2. Let (G, I, O) be a geometry. A collection C of (possibly trivial) directed paths in G is a path 
cover of (G, /, O) if 

(«). each V 6 V{G) is contained in exactly one path (i.e. the paths cover G and are vertex-disjoint); 

(a), each path in C is either disjoint from / , or intersects / only at its initial point; 

(Hi), each path in C intersects O only at its final point. 

The successor function of a path cover C is the unique / : O*^ — > I'^ such that y = f{x) if and only li x ^ y 
is an arc of C . If a function / : O"^ — > P is a successor function of some path-cover of (G, 7, 0) , we call / a 

successor function of (G, /, O). 
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If a geometry {G,I,0) has a causal flow (/, =4) , the maximal orbits of the successor function / define a path 
cover for (G, /, O) , which allow us to consider the causal flow in terms of vertex-disjoint paths in G : 

Theorem 3 [7, Lemma 3]. Let (/, be a causal flow on a geometry {G,I,0) . Then there is a path cover 
Vf of (G, I, O) whose successor function is f . 

Given that the successor function of a causal flow for (G, /, O) induces a path cover, one might think of also 
trying to obtain a causal flow from the successor function of a path cover. There is an obvious choice of binary 
relation for a successor function / : 

Definition 4. Let / be a successor function for (G, /, O) . The natural pre-order^ =4 for f is the transitive 
closure on V{G) of the conditions 

X =4x ; X =4 f{x) ; y ~ f{x) =^ x =4 y , (2) 

for alla;,y € V{G) . 

If =^ is a partial order, it will be the coarsest partial order such that (/, =4) is a causal flow. However, it is easy to 
construct geometries where =^ is not a partial order. Figure 1 illustrates one example. For any choice of successor 
function / on this geometry, (Fm) forces either ao =4 o-i =4 0.2 =4 «o or ao ^ ai ^ a2 ^ ao to hold. Because oq , 
ai , and 02 are distinct, such a relation =<; is not antisymmetric, so it isn't a partial order. 



ao 



ai 



02 




ao =^ bo 

ai =4 bi 

0-2 ^ b2 

ao ^ ai =4 a2 ^ ao 



I O 

Figure 1: A geometry with a successor function / : O'^ — > /"^ , but no causal flow. 



In the example above, we have a cycle of relationships induced by condition (Fm). The following definitions 
characterize when such cycles of relationships occur. 

Definition 5. Let (G, /, O) be a geometry, and a family of directed paths in G . A walk W = uoUi ■ ■ - u^ 
is an influencing walk^ for if it is a concatenation of zero or more paths (called segments of the influencing 
walk) of the following two types: 

• xy , where a; — > y is an arc of .F ; 

• xzy , where a; ^ z is an arc of J- and yz G E{G) . 

A vicious circuit for ^ is a closed influencing walk for with at least one segment. 

Theorem 6 [7, Lemma 9]. Let C be a path cover for {G,I,0) with successor function f , and let =4 be the 
natural pre- order of f . Then x ^ y if and only if there is an influencing walk for C from x to y . 



Given that we want to forbid cycles of relationships for the natural pre-order =<; , we are then interested in the 
following restriction of path covers: 



•^A prc-ordcr is a binary relation which is reflexive and transitive, but not necessarily antisymmetric. 
■^These are closely related to walks which alternate with respect to T: see Section 3.1.1. 
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Definition 7. A path cover C for (G, I, O) is a causal path cover if C does not have any vicious circuits in G . 

Theorem 8 [7, Theorem 10]. Let {G,I,0) be a geometry with path cover C , f be the successor function of 
C , and ^ he the natural pre- order for f . Then C is a causal path cover if and only if ^ is a partial order, which 
occurs if and only if (/, =4) is a causal flow for (G, /, O) . 

By characterizing causal flows in terms of causal path covers, we can make use of the following result: 

Theorem 9 [7, Theorem 11]. Let {G,I,0) be a geometry such that \I\ = \0\ , and let C he a path cover for 
{G,I,0) . If C is a causal path cover, then C is the only maximum collection of vertex-disjoint I - O dipaths. 

Then, if |/| = \0\ and {G,I,0) has a causal flow, there is a unique maximum-size collection of vertex-disjoint 
I - O paths, and that collection is a causal path cover which allows one to reconstruct a causal flow. Taking 
the contrapositivc, if wc can find a maximum-size collection of vertex-disjoint paths from 7 to O which is not a 
causal path cover, then (G, /, O) does not have a causal flow. 

3 An efficient algorithm for finding a causal fiow when |/| = \0\ 

Using Theorems 8 and 9 when |7| = \0\ , we can reduce the problem of finding a causal flow to flnding a 

maximum-size family of vertex-disjoint I O paths in G . Given such a family of paths J-' , wc may then verify 
that the resulting family forms a path cover for G , obtain the successor function f oi !F , and attempt to build 
a causal order compatible with / . We illustrate how this may efliciently be done in this section. 

Implementation details. For the purpose of run-time analysis, I fix here conventions for the data structures 
used to implement graphs, paths, and sets throughout the following algorithms. 

• We will assume an implementation of graphs and digraphs using adjacency lists for each vertex x (in the case 
of digraphs, using two separate lists for the arcs entering x and those leaving x). Such an implementation 
can be easily performed in space 0{m) , where m is the number of arcs/edges, assuming a connected 
(di-)graph.^ 

• Sets of vertices are considered to be implemented via arrays storing the characteristic function of the set. 

Wc may assume without loss of generality that these are also used to perform bounds-checking on arrays 
which are used to implement partial functions on V{G) , such as successor functions / : O'^ — > I'^ . 

• Collections of vertex-disjoint di-paths in a graph G will be implemented as a set V{J^) indicating for each 
X e V{G) whether x is covered by , and a digraph containing all of the arcs ofj^. As well, functions prev 
and next will be defined for all vertices in /'^ (respectively, O*^) covered by which returns the predecessor 

(respectively, successor) of a vertex covered by T . 

Throughout some of the algorithms below, a family of vertex-disjoint paths may be transformed into to a 
graph where a single vertex has out-degree 2 , but every other vertex has out-degree at most 1 , and every 
vertex has in-degree at most 1 . So long as these bounds are maintained, determining whether a vertex is 
covered by J-", whether an arc is in T , and adding/deleting arcs from can be done in constant time. 
As well, the function prev will be well-defined so long as the in-degree of the graph representation of is 
bounded by 1 . 

''This holds, in particulcir, for graphs corresponding to one-way patterns implementing unitary operations which are not tensor- 
product decomposable. 
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3.1 Efficiently finding a path cover for {G,I,0) 



Given a geometry (G, /, O) , we are interested in obtaining a maximum-size family J- of disjoint I-O paths in G 
in order to test whether it is a causal path cover. This is known to be efficiently solvable. 

Problems involving constructing collections of paths with some extremal property in graphs are usually solved 
by reducing the problem to a problems of network flows on digraphs: algorithms for such problems have been 
very well studied. (Section 4.1 of [7] outlines an algorithm of this kind to find a maximum-size family of disjoint 
1-0 paths.) However, in order to present a solution which does not assume any background in graph-theoretic 
algorithms, and also in order to reduce the number of auxiliary concepts involved in the solution, I will present 
an algorithm not explicitly based on network flows. A dividend of such a presentation is that it highlights the 
relationship between influencing walks and walks which alternate with respect to a collection of disjoint paths, 
which was alluded to in Definition 5. 

3.1.1 Alternating and augmenting walks 

Definition 10. Let I,0 C V{G) . A collection of vertex-disjoint paths from / to O is proper if its' paths 
intersect / and O only at their endpoints. 

A collection of k vertex-disjoint I-O paths of is necessarily proper when |/| — \0\ — k . We would like to arrive 
at such a maximum-size collection by producing successively larger proper collections of vertex-disjoint paths. 
To so so, we will use results of graph theory pertaining to Menger's Theorem. The basic approach present is 
outlined in Section 3.3 of [9]. 

Definition 11. For a family of vertex-disjoint directed paths from / to O , a walk W = uqUi ■ ■ ■ U( in G is 
said to be pre- alternating with respect to T if the following hold for all < j, fc ^ £ : 

(i). T does not contain Uj Uj+i as an arc; 
(ii). if Uj — Uk and j ^ k , then Uj is covered by J- ] 

{Hi), if Uj is covered by J- , then either Uj — > Uj-i or Uj4_i — > Uj is an arc of T . 

W is said to be alternating with respect to J- if W is pre-alternating with respect to J-' , and uq is an element of 
/ not covered hy J- . is an augmenting walk for T ilW alternates with respect to J- , and ui £ O . 

Figure 2 illustrates two pre-alternating walks for a family J- of vertex-disjoint paths in a geometry (G, /, O) . 



y-O M4 




Figure 2: Two examples of a walk W (hollow arrows) which is pre-alternating with respect to a collection of vertex-disjoint paths 
from I to O (solid arrows). In both examples, circled vertices are entry points of W into (see Definition 13). 

The relationship between influencing walks and pre-alternating walks is most clear for a path cover C of (G, /, O) , 
in which case an influencing walk for C is the reverse of a walk which is pre-alternating for C . As we will see 
in the next few pages, pre-alternating walks describe ways in which different families of disjoint paths from / to 
O are related to each other: this is essentially the reason why a vicious circuit (i.e. a closed influencing walk) 
exists for a path cover whenever there is a second family of disjoint I-O paths of the same size. 

^The solution presented here can be easily related to the solution via network flows, but a small amount of additional work must 
be done in order to stay in the context of collections of disjoint paths, rather than disjoint paths, cycles, and walks of length 2 . 
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First, we will show that augmenting walks for T are always present if \T\ < \I\ — \0\ , and if there is a family 
of vertex-disjoint I - O paths of size k : 



Theorem 12. Let G be a graph, and 1,0 C V{G) with \I\ = \0\ = k . Let T he a collection of vertex- disjoint 
I - O paths with \J-\ < k , and i G L be a vertex not covered by T . Lf there is a collection C of vertex-disjoint 
dipaths from L to O with \C\ = k , then there is an augmenting walk W for T starting at i which traverses each 
edge of G at most once, and where i is the only input vertex in W not covered by J- . 

Proof — Suppose G contains a collection C of fc vertex-disjoint I - O dipaths, let be some proper collection 
of vertex-disjoint I - O dipaths of size less than k , and let /' 7^ be the set of input vertices not covered by . 
Let us say that a vertex v G V{G) is an incidence point of C and if w is covered by both C and J- , and there 
is a vertex w which is adjacent to u in a path of C but which is not adjacent to v in any path of J- . Let 3 be the 
set of incidence points of C and T , and let S be a di-graph with V{9) = /' U 3 U O , and {x ^ y) E ^(S) for 
x,y G V{Q) if one of the following applies: 

• there exists a vertex z G such that 

(z). x and z lie on a common path P in C , where z is the next incidence point in P after x , and 
(ii). y and z lie on a common path P' in J- , where z is the next incidence point in P' after y ; 

• X and y lie on a common path P in C , there are no incidence points on P after x , and y G O . 




Figure 3: Two families of vertex-disjoint I - O paths in a 
graph: one family C with k paths (hollow arrows), and one 
family with < k paths (solid arrows) . Circled vertices are 
the incidence points of C and . Dashed lines are the other 
edges of the graph. 



Figure 4: The digraph S obtained by applying the construc- 
tion above to Figure 3. Dashed arrows represent the edges 
from path segments belonging to either C or in the origi- 
nal graph; thick black arrows are the actual arcs of S , which 
are induced by those path segments. 



Because both C and J-' are vertex-disjoint collections of paths, it is easy to show that the maximum in-degree 
and out-degree of S are both 1 . Thus, S consists of vertex-disjoint di-cycles, walks of length 2 , isolated vertices, 
and directed paths. 

Because each u g /' is not covered by a path of J- , and is not preceded by any vertices in it's respective path of 
C , it has in-degree in S . Then, each element of /' is at the beginning of a maximal dipath in S . Furthermore, 
each vertex in w G /' U 3 has out-degree 1 : if P G C is the path covering v , either there are no incidence vertices 
after v on P , in which case there is an arc u — > y for the vertex y G O at the end of P ; or if we let z G 3 be 
the first incidence vertex following w on P , we will have z ^ I , in which case there will be an incidence vertex 
w which precedes z on some path of JF, because all input vertices covered by T are incidence points. Thus, any 
maximal di-path in 9 must end in O . Then, for each i £ L' , there is a dipath from i to some element of O in 
the graph S . 

Consider any vertex i G I' , and let uq ui —>■■■■ ^ ue be the dipath in 9 from i to O . Let P G C and P' E J- 
be the paths containing ui-i : from {ui-2 G ^(9) , we know that there is an incidence vertex after itf_i 

in the path P' . Note that ui is either not covered by any path oi !F , or it occurs at the end of a path of !F and 
is not followed by any vertices on that path; then, (u£_i U() G ^(9) implies that there are no incidence points 
on P after Uf_i . Then, the arcs leaving in P and P' must be different: the fact that no incidence point 
follows Uf-i in P then implies that no path of J- intersects P after ui_i . In particular, ui is not covered by J- . 
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Because uq G /' and G O arc both not covered by , we may construct an augmenting walk W for in the 
original graph G , as follows. If ^ = , we let be the trivial path on uq , which is an augmenting walk for . 
Otherwise: 

• For each j G — 1] , let vj be the next incidence point after Uj on the path Pj G C containing Uj . (This 
Vj will then also be the next incidence point after u^i on the path Pj G containing Uj^i .) 

• Let Pj be the segment of Pj from Uj to Vj , and be the reverse of the segment of Pj from Uj^i to Vj . 

• Finally, let P^-i be the path segment in C from U£-i to ue . 

Then, define W = uoPoVoPqUiPi ■ ■ ■ ue-iPe-iUg . We may show that W is an augmenting walk for : 

(i). Each path Pj is internally disjoint from , because they are sub-paths of elements of C, and do not 
contain any incidence points in their interiors. Then, none of the arcs of Pj arc arcs of !F for any j G [£] . 
Also, all of the arcs of the paths Fj are the reverse of arcs of : they do not contain arcs of ^ either. 
Then, none of the arcs of W are arcs of . 

(m). Because uq ^ • ■ ■ ^ ue is a directed path in S , we have Uj ^ Uk ■ Because each path Pj and P- is 
uniquely determined by Uj for j G [£—1] , those sequences of vertices can also occur only once each. Each 
interior vertex of Pj or Pj can only occur in a single path of C or .7-" , between two consecutive elements 
of 7' U J U O on that path : then, because each segment Pj and Pj only occur once in W , each interior 
vertex of those segments also occurs only once in W . 

Thus, if any vertex x occurs more than once in W , x miist be an end point of some path Pj or Pj . 
Aside from Pq and Pe-i , both end-points of each such segment has in-degree 1 and out-degree 1 , so 
they cannot be elements of either I' or O . Then, any vertex which occurs more than once in W is an 
clement of 3 , and is therefore covered by . 

{in). The only points in W which are covered by are the vertices of the paths Pj for j G — 1] , which are 
all at the beginning or the end of arcs in W which are the reverse of arcs of . 

Thus, W is an augmenting walk for T . Furthermore, because each edge of G is contained in at most one segment 
Pj or Pj , each edge occurs at most once in W . Finally, because elements of I' have in-degree in S and do 
not occur in the segments Pj or Pj , any input vertices other than i = uq which occur on W must be incidence 
points, which means they are covered by J-" . Thus, there is a proper augmenting path for of the desired type 
starting at i G / . □ 

The above Theorem illustrates how we can build an augmenting walk for from a collection of disjoint I - O 
paths which covers / and O . If we impose restrictions on the type of augmenting walk we consider, we may also 
efficiently do the reverse. The restriction we are interested in is the following: 

Definition 13. Let W = uqUi ■ ■ - uehe a walk which which is pre-alternating with respect to . 

• An entry point of W into ^ is a vertex Uj which is covered by , where either j = or Uj — > is not 

an arc of J- . 

• The walk W is monotonic if, for every path P & T and for any indices ^ h < j < i such that Uh and Uj 
are both entry points for W into which lie on P , Uh is closer to the initial point of P than than Uj is. 

• W is a proper pre-alternating walk if W traverses each edge at most once, each input vertex in W (except 
possibly Uq) is covered by ^, and W is monotonic. 

We will be most interested in proper augmenting walks, which are useful in increasing the size of proper collections 
of 7 - O paths. The sort of augmenting walk that is guaranteed by Theorem 12 is almost a proper augmenting 
walk, and merely lacks a guarantee of monotonicity. However, the following Lemma shows that we lose no 
generality in imposing monotonicity as a condition: 
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Lemma 14. Let G be a graph, and I,0 ^ ■ Let J- be a collection of vertex-disjoint L - O paths, and let 

W be an augmenting walk for T from i ^ I to ui & O . Then there is a monotonic augmenting walk W for T 
from i to LO . 

Proof — Let W be given by W — uq ■ ■ ■ ui , where uq ~ i and ui ~ lu . For any path P E J- , and two 
entry points Uh and Uj of W into T , let us say that {uh,Uj) is a reversed pair if h < j but Uj is closer to the 
initial point of P than Uh ■ We will produce a monotonic augmenting walk by recursively reducing the number 
of reversed pairs of W . 

• If has no reversed pairs, then W is already monotonic, in which case we may let W = . 

• Suppose that (uh, Uj) is a reversed pair of W . Then h < j , but Uj is closer than to the initial point of the 
path Q £ J- which covers both of them. Note that Ui — lo is not covered by J- , and so is not on the path Q : 
then, let j' G [£] be the smallest index such that Wj'+i is not on Q . Let Q = qo ■ ■ ■ QaQa+i • ■ • Qb-ilb • • ■ , 
where qa — Uji and qf, — . Then, let 

W' ^Uo-- ■ Uh-iqbqb-iqb-2 ■ ■ ■ qa+iqaUj' + l ■■ - Ui . 

From the fact that W is an augmenting walk for it is easy to show that W' is also an augmenting walk 
for . As well, the entry points of W into J- are a subset of the entry points of W into !F , in which case 
the reversed pairs of W are also a subset of the reversed pairs of W ; and W' does not have {uh, uj) as a 
reversed pair. Then, W' has strictly fewer reversed pairs than W . 

Because is a finite walk, it can have only finitely many reversed pairs; then, by recursion, we may construct 
a monotonic augmenting walk W for J- from i to u . □ 

Corollary 15. Suppose \L\ = \0\ — k , !F a proper collection of vertex- disjoint I - O paths in G with \ J- \ < k , 
and let i £ I be a vertex not covered by T . Lf there is a collection C of vertex-disjoint dipaths from L to O with 
\C\ = k , then there is a proper augmenting walk W for J- starting at i . 

Proof — Theorem 12 and Lemma 14. □ 

For proper augmenting walks, the reason for requiring that no edge is traversed twice is essentially to help 
construct efficient algorithms for finding them, which we consider later. The requirements that the only input 
vertex in the walk which is not covered by T , and that it be monotonic, are essentially chosen to allow us to use 
augmenting walks to increase the size of a proper collection of vertex-disjoint paths to cover exactly one more 
input vertex. We may do this using the following operation: 

Definition 16. Let be a proper collection of vertex-disjoint I - O dipaths in G, and W he & proper 
augmenting walk for T . Then, !F ®W denotes the collection of directed paths which are formed by those arcs 
X y which belong either to or a path of J- , and for which y —^ x is not an arc of either W ov J- . 
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The collection J- (SW described above is formed by the usual procedure for augmenting a network-flow with an 
augmenting path: one can think of forming J^(BW hy "adding" together the arcs of !F and W , and "cancelling" 
them whenever they point in opposite directions on a single edge. 

Lemma 17. Let J- be a proper collection of vertex-disjoint I - O dipaths in G , and W he a proper augmenting 
walk for T . Then T ®W is a proper collection of vertex-disjoint I - O dipaths, with \ ®W\ = \ J^\ + 1 ; and 
the input vertices covered by ®W are those covered by T and W together. 

Proof — We induct on the number of times r that the walk W intersects the paths oi T . If r = , then 
J- ®W = T U {VK}, and the inputs covered by ^ © W are clearly those covered by ^ or by . Otherwise, 
suppose that the proposition holds for all cases where the augmenting walk intersects the paths of its' respective 
collection fewer than r times. 

Lot W be given by W = U[)U\ ■ ■ ■ UaUa+x ■ ■ ■ Uf^-iUf, ■ ■ ■ ui , where none of the points Uj are covered by T for j G [a] , 
and where uj^i — *■ Uj is an arc of J' for all a ^ j < b . Let Q G he the path containing Ua through Uf, : in 
particular, let Q = qoqi ■ ■ ■ gc^ofi • • • 9d-i9d ■ ■ - Qm, where Qc = Ub and qa = Ua ■ Then, we may define 

Q' = UqUi ■ ■ ■ UaQdH ■■■qm , W' = q^qi ■ ■ ■ qcUhi-i ■■■ui : 

then Q' £ (BW , and W is an augmenting walk for .F' = \ Q) U {Q'} which intersects the paths of J^' 
fewer than r times. Because is a proper collection of vertex-disjoint I - O paths, Q' only intersects / and 
O at its' endpoints, and Q' does not intersect any paths of \ Q , T' is proper. Similarly, because Q only 
intersects / at q^ and because W only intersects / at uq and at input vertices covered by T , W does not cover 
any inputs except those covered by T . Because W doesn't traverse any edges twice, and all of the other entry 
points qh of W into T on the path Q have ft. > c by the monotonicity of W , W itself does not traverse any 
edge twice. Finally, all of the entry points of W into T are also entry points of W into T' , except for Ua : all 
the other are left unaffected, including the order in which they occur. Then W' is monotonic, so that W is a 
proper augmenting walk for T' . 

By the induction hypothesis, T' ®W is a proper collection of vertex-disjoint paths from J to O , with \T' © W'\ = 
+ 1 = I J^l + 1 . Also by induction, the input vertices covered by T' © W are those covered by T' or by 
W . Because W covers the input go , and T' covers all inputs covered by or by except for go , J^' ® W 
then covers all vertices covered by .7-" or by . Finally, note that the set of arcs from and W together only 
differs from the set of arcs from T and W together by the absence of the arcs Uj Uj+i from W and the arcs 
Uj+i ^ , for a < j < 6 , which oppose each other. We then have J^' © W = J^®W : thus, |.F ffi PF| = -|- 1 , 
and ®W covers the input vertices covered either by or by W . □ 

3.1.2 An efficient algorithm for finding a proper augmenting walk 

Algorithm 1 determines if a vertex supports a suitable proper pre-alternating walk W with respect to T , and 
compute J-(SW if one is found. Using it, we may find proper augmenting walks for J- by performing a depth-first 
search along proper alternating walks W for J- in an attempt to find one which ends in O . 

Theorem 18. Let {G,L,0) he a geometry with \I\ = \0\ = k . T a proper collection of fewer than k vertex- 
disjoint paths from I to O , iter a positive integer, i € I a vertex not covered by J- , and visited : V{G) — > N 
with visited(a;) < iter for all x e V{G) . Then Augment Search halts on input (G, 7,0, .F, iter, visited, i) . 
Furthermore, let {T, visited, status) = AugmentSearch(G, 7, 0,.F, iter, visited, i) . 

(i). If status = fail , then there are no proper augmenting walks for T starting at i ; 

(n). // status = success , then T is a proper family of vertex-disjoint I O paths of size \J^\ + 1 which 
covers i and all input vertices covered by T , and visited(a;) ^ iter for all x G V(G) . 

Proof — Let G, 7, O, T ^ and iter be fixed as above. Throughout the proof, we will consider chains of 
recursive calls to AugmentSearch. One invocation of Augment Sear oh. is the daughter of a second invocation if 
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Algorithm 1: AugmentSearch(G, /, O, JT, iter, visited, — searches for an output vertex along pre- 
alternating walks for J- starting at v , subject to limitations on tlie end-points of tlie search paths. 
Require: {G, I, O) is a geometry. 

Require: is a specification for a vertex-disjoint family of I-O paths. 
Require: iter is a positive integer. 
Require: visited is an array V{G) — > N . 
Require: v G V{G) . 



1: visited(w) ^ iter; 

2: ii V E O then return(jF, visited, success). 

3: if V € V{J^) and v ^ I and visited(prev(jF, v)) < iter then 

4: (.F, visited, status) <— AugmentSearch(G, I, O, iter, visited, prev(.F, u)) ; 

5: if status = success then 

6: .F <— RemoveArc(^, prev(.7-', v) ^ v); 

7: return (.F, visited, success) . 

8: end if 

9: end if 

10: for all w ~ ?; do 

11: if visited('u;) < iter and w ^ I and — > w) ^ ^(■^) then 
12: if w ^ V{T) then 

13: (.F, visited, status) ^ AugmentSearch(G, /, O, .F, iter, visited, w) ; 

14: if Status = success then 

15: •!— AddArc(.F', u — » w) ; 

16: return (^, visited, success) . 

17: end if 

18: else if visited(prev(jF, w)) < iter then 

19: (.F, visited, status) ■s^ AugmentSearch(G, /, O, J-, iter, visited, prev(.F, w)) ; 

20: if status = success then 

21: RemoveArc(J-", prev(^, w) — > w) ; 

22: JF <— AddArc(.F', v ^ w)\ 

23: return (.F", visited, success) . 

24: end if 

25: end if 

26: end if 

27: end for 

28; return (^, visited, fail) . 



the first invocation was performed as a step of the second invocation; if one invocation is related to a second 
invocation by a sequence of daughter-relationships, we will call the second invocation a descendant of the first. 

At any stage in a particular invocation of AugmentSearch, we will refer to the ordered pair (visited, v) as the 
data pair of the invocation, where v is the final parameter of the input, and visited the second last parameter, 
including any changes which have been made to it during the invocation. (Though the input parameters of 
AugmentSearch include G, /, O, and iter, we will occasionally refer to data pairs as the input of an 
invocation of AugmentSearch.) When an invocation of AugmentSearch has a data pair (visited, v) and makes a 
daughter invocation, we may describe that invocation as being "daughter invocation for (visited, f)" ; similarly, 
a daughter invocation for (visited, w) or the descendant of one is a "descendant invocation for (visited, v)". 

We define a probe walk W for an ordered pair (visited, i;) to be a proper pre-altcrnating walk starting at v 
such that, for all vertices x in the walk, visited(x) — iter only if x is at the beginning of "W and x & O only 
if x is covered by .F or a; is at the end of W . Then, we let i?(visited, v) be the set of vertices x G V{G) which 
end-points of probe walks for (visited, v) . We will reduce the problem of determining whether there is a proper 
augmenting path for .F passing through w to a question of the existence of whether there is an output vertex in 
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i?(visited, v) , for visited restricted in a manner described below. 

A canonical walk for a data pair (visited, u) is a proper alternating walk W with respect to , such that the 
following all hold: 

(i). visited(a;) ^ iter for all x e V{G) . 
{ii). V is the end-point of W . 

{Hi), for all vertices x on W , if visited(a;) < iter , then either x = v , or x occurs exactly once in W and is 

an entry point of W into J- . 

{iv). for any path P of J- , and x G V{P) which is not in W , visited(x) < iter if and only if either (a) there 
is no entry point of W after x on the path P , or (b) there is exactly one entry point pofW after x on 
the path P , and v lies on P strictly between x and p . 

A data pair (visited, t;) is itself canonical if it has a canonical walk. We will be interested in the behaviour 
of AugmentSearch on canonical inputs. (Note that the input described in the statement of the Theorem is a 
special case.) We will show that AugmentSearch essentially performs a depth-first traversal of i?(visited, 
along probe walks for (visited, v) in an attempt to find an output vertex. If it succeeds, it has traversed a 
proper augmenting walk W for .F , and can construct ® W . 

Suppose is a canonical walk for a data pair (visited, v) . It is easy to show that if we extend W to a. longer 
walk W' = Wv'WiW2 ■ ■ -wpf for some N ^ 1, and W is a canonical walk for some data pair (visited*, Wjv) ; 
then is a canonical walk for (visited*, u) . We will use this fact frequently in the two Lemmas below. 

Lemma 18-1. Suppose that (visited, t;) has a canonical walk W . If R{visited,v) does not contain any 
output vertices, AugmentSearch halts on input da,ta (visited, w) , with output value {T, visited, fail) ; where 
visited differs from visited only in that visited(a;) = iter for all x S i?(visited, v) , and where (visited, v) 
also has the canonical walk W . 

Proof — We will proceed by induction on the length i of the longest probe walk for (visited, v) . Regardless 
of the value of £, line 1 transforms the data pair (visited.?;) to (visited^^^ ?;) . where visited^^^ differs from 
visited in that visited^^''(i;) = iter; then, any canonical walk for (visited,!;) is also a canonical walk for 
(visited^^^t;) . As well, it cannot be that v £ O: then the condition on line 2 will not be satisfied. 

If £ = , the condition on lines 3 cannot be satisfied, and the condition of line 11 is not satisfied by any neighbor 
w ^ V. Then, Hnc 28 will ultimately be executed, returning (^, visited^^\ f ail) . Because i?(visited, u) = 

{v} , the proposition holds in this case. 

Otherwise, suppose i > 0, and that the proposition holds for canonical data pairs whose probe walks all have 
length less than i . Consider the vertices which may be the subject of a daughter invocation of AugmentSearch: 

1. li V ^ I and v is covered by , and z is the predecessor of v in the paths of , then (visited,!;) has 
probe walks starting with the arc u — > 2; if and only if visited(^^(2;) = visited(2:) < iter. If this holds, 
then a daughter invocation of AugmentSearch with input data (visited^^^ z) is performed. 

In this case, note that (visited*^^\ 2;) has probe walks ending in O only if (visited,!;) does; then 
i?(visited^^\ 2) is disjoint from O , and all of the probe walks of (visited^^\ 2;) are strictly shorter than 
those of (visited, w) . Let W''^^ = Wvz: because W is a canonical walk, z G l^(M^) only if z occurs only 
once in W and is an entry point of W into , in which case the edge vz is never traversed by W . Then, it is 
easy to show that W^^^ is a canonical walk for (visited^^-* , z) . By the inductive hypothesis, AugmentSearch 
will halt on input (visited^'^-' , z) and return a value (JF, visited^^' , f ail) , where visited^^-* differs from 
visited^^) only in that visited(^^(a:) = iter for all x € i?(visited(^\ ^) C R(yisited,v) , and where 
W^^^ is a canonical walk for (visited^^^ 2;) . Then, is a canonical walk for (visited^^^, v) . 

Otherwise, if visited^^^ (z) = iter , if i; G 7 , or if F" does not cover v , let visited^^^ = visited^^^ ; W is 
a canonical walk for (visited^^^, u) in this case as well. 

2. Suppose that at some iteration of the for loop starting at line 10, the data of AugmentSearch is a data 
pair (visited^''^ v) for which 14^ is a canonical walk, visited('')(!;) = iter, and w is a neighbor of v 
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satisfying the conditions of lines 11 and 12. Note that (visited^''^ w) has probe walks ending in O only if 
(visited, v) does; then, i?(visited('') , w) is disjoint from O , and all of the probe walks of the former are 
strictly shorter than those of the latter. Let W^'^^ = Wvw ; because w is not covered by T , the fact that 
visited*^'*'' (w;) = visited(ij;) < iter implies that w does not occur in W . Because w ^ I, we know that 
W^^'^ is a proper alternating walk. In particular, M^C*) is a canonical walk for (visited^''-', w) . 

By the inductive hypothesis, Augment Search will then halt on input (visited^''^ w;) and return a value 

{J-, visited('*+-'^^ , fail) , where visited'^''+-'^^ differs from visited^'') only in that visited('*+-'^^ (.r) = iter 
for all X G i?(visited('*\ w) C i?(visited, ti) , and where W^'^^ is a canonical walk for (visited*^''+^\ 
Then, is a canonical walk for (visited(''+^\t;) . 

3. Suppose that at some iteration of the for loop starting at line 10, the data of Augment Search is a data pair 
(visited'^''^ u) for which is a canonical walk, visited^''^ (u) = iter , and w is a neighbor of v satisfying 
the conditions of lines 11 and 18. Then, w is covered by a path P and has a well-defined predecessor 
^; in P . Let W^C*) = Wvwz : this is an alternating walk with respect to . 

The walk M^C*) is monotonic only if w is further from the initial point of the path P G than any entry 
point of on P . If P contains no entry points of W into , this is satisfied. Otherwise, let y be the final 
entry point of W into P . 

• Suppose that v is not covered by P . Because (visited^'^^ w) is canonical, every vertex x on the path 
P from the initial point up to (but possibly not including) y has visited(''^(a;) = iter. Because 
visited^''' (2;) < iter, z is at least as far along P as y is; then, w is strictly further. Thus, W^'^'> is 
monotonic. 

• If w is covered by P, then every vertex x on P with visited(x) < iter either is at least as far as 
y on P , or has the property that y is the only entry point between x and the end of P , and that 
V lies between x and y . However, if there are more than zero vertices of the second type, then v 
has a predecessor z in P with visited(2) < iter. Then, from the analysis of part 1 above, all 
vertices x which precede v in P with visited(a;) < iter are in ii(visited(^^, 2;) , and thus have 
visited^''^ (.t) = visited^^^ (x) = iter. Then, T^^C*) is monotonic if and only if w is further along P 
than y , which reduces to the analysis of the preceding case. 

Because visited('')(2;) < iter , either z is not inW , ot it occurs exactly once as an entry point of W into 
. Because visited^''^(w) < iter and w is further along P than any entry point ofW , w does not occur 
in W at all. Then, neither vw nor wz are traversed by , in which case W^^^ is a proper alternating walk. 

In particular, it is a canonical walk for (visited^'*^ z) . 

Again, (visited^'*\ z) has probe walks ending in O only if (visited, z;) does; then, P(visited^'*\ z) is 
disjoint from O , and all of the probe walks of the former are strictly shorter than those of the latter. 
By the inductive hypothesis, AugmentSearch will then halt on input (visited^'*^ z) and return a value 
(^, visited'-'*"'"^-', f ail) , where visited^'*^^-' differs from visited'-'*^ only in that visited'-'*"'"^-' (x) = iter 
for all X G P(visited^''\ 2:) C P(visited, v) , and where W^'^^ is a canonical walk for (visited^''+^^, 0) . 
Then, is a canonical walk for (visited('*+^),t;). 

By induction on the number of neighbors w ~ v satisfying the conditions of lines 11, 12, and 18, the data 
(visited, t;) when the for loop terminates and line 28 is executed will be a nearly canonical pair, and visited 
differs from visited only on elements of P(visited, v) . 

It remains to show that visited(a;) = iter for all x G P(visited, v) . We have shown this already foY x = v; 
then, let r G P(visited, t;) \ {v} . By definition there is a probe walk W for (visited, u) ending in r. The 
vertex w immediately following w on W will be either tested on line 3 or line 11 as a neighbor of v ; then, there 
exists indices h' such that W is not a probe walk of (visited*^'' \v) . Let h > Qhe the largest integer such that W 
is a probe walk for (visited'''^ u) : then, there are vertices x ^ v inW such that visited(''+^^(a;) = iter . Let 
y G P be the last such vertex in W , let W be the segment of W from y onwards: then r G P(visited(^+-'^\ y) . 
Because visiteS'^^^\y) = iter, there must have been a descendant invocation for (visited' ''^ u) which had 
input data (visited*, y) for some function visited* : it is not difficult to show that visited* (y) < iter. By 
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induction on daughter invocations using the analysis above, we may show that (visited*, y) is a canonical data 
pair with probe walks strictly shorter than £ : then, for all x e i?(visited*, y) , we have visited'-''+^^(a;) — iter . 
However, because visited(''+^^ (x) < iter impHes visited* (x) < iter , and because all vertices x £ F(W) after 
y have visited(''+^'(a;) < iter, W is a probe walk for (visited*,?/) . Then, we have visited(''+^) (r) = iter. 

Because visited(a;) = iter if and only if visited('*)(a;) = iter for some /i ^ 1 , we then have visited(r) = iter 
for any r G i?(visited, v) . By induction, the Lemma then follows. 

Lemma 18-2. Suppose that (visited, u) is a canonical data pair. // i?(visited, u) contains an output vertex, 
AugmentSearch halts on input data (visited,?;) , with output value (.F © W, visited, success) ; where W is a 
probe walk for (visited, i>) ending in O , and visited differs from visited only in that visited(a;) = iter 
only for x in some subset of R{v±s±ted,v) . 

Proof — We induct on the length f e N of the longest probe walk for (visited.?;) ending in O . If £ = , 
then V G O , and the result holds trivially. Otherwise, suppose £ > Q and that the result holds for those canonical 
data pairs (visited*, x) which have probe walks of length less than i ending in O . 

Let be a canonical walk for (visited, ?;) . Consider the sequence of vertices wi ,W2 ,■ ' ' > wm which are tested 
(either on line 4, line 13, or line 19) in the course of the invocation of AugmentSearch. We let visited*^^) differ 
from visited in that visited'^^-' (?;) = iter, and from this define visited^-' •* for j > 1 by letting visited*^-'+^-' 
be the second component of the output of the daughter invocation with input data (visited^-'^Wj) . (If the 
daughter invocation with data pair (visited^^^ ?«m) halts, this sequence extends to visited^^+^^ .) 

Let W be a probe walk for (visited, ?;) which ends at a vertex uj <= O , and let 1 ^ ^ + 1 be the largest 
integer such that ii(visited(''^, ??;j) is disjoint from O for all j < . If = M + 1 , this means that the 
invocation of AugmentSearch on input data (visited^^^ ?i;M) halted with fail in the final part of its' output 
value, and that there are no neighbors of ?« ?; which can satisfy the conditions of lines 11, 12, and 18 (due 
to the choice of M as the length of the sequence of daughter- invocations). However, we may show by induction 
that for all 1 < j < A'' , W is a probe walk for (visited^-'^ ?;) , which is canonical: 

• This follows immediately for j = I , because visited*^^) only differs from visited at v , and thus has W 
as a probe walk and as a canonical walk. 

• Suppose for some 1 < j < M that W is a canonical walk for (visited^-'^t;) , that W is a probe walk 

for (visited^-'^ ?;) , and that visited^-'^ (?;) = iter. Then we can extend W to a. canonical walk W^^^ 
for (visited^-'^ lOj) : either by setting W^^^ = Wvwj in the case that {wj —>■?;) S A{J^) or Wj is not 
covered by J^, or by setting W^^^ = Wvzwj where {z wj) e A{T) otherwise. Because R{visited^^\ Wj) 
contains no output vertices, by Lemma 18-1 we know that visited^^"*"-*^^ differs from visited^-'^ only on 
i?(visited*^-^\ ??;j) and that W^^^ is a canonical walk for (visited(''+^\ ?i;j) . Then, is a canonical walk 
for (visited^-'"'"^^ ?;) . 

For any vertex a; in W, the sub-path Wx from x to w is a probe walk for (visited^-'^x) of length less 
than i. If W has a non-trivial intersection with i?(visited(^\ ??;j) , then some vertex x G T^(W) is the 
first such vertex which is given as part of an input data pair (visited*, x) for a descendant invocation 
for (visited^-'^, Wj) . By induction on the recursion depth from ?; to a; , we may show that there is then a 
probe walk W* for (visited'^-'),?;) ending at x, and that VFW* is a canonical walk for (visited*, x) ; and 
precisely because x is the first vertex of W which is visited in a descendant invocation for (visited*^-*-*, Wj) , 
we know that visited* (y) < iter for all y G VCW) \ {?;, x} . Then, VJ^ is a probe walk for (visited*, x) , 
and by the induction hypothesis, this invocation of AugmentSearch then terminates with success as the 
last part of its' output. Again by induction on the recursion depth, we may also show that the invocation 
of AugmentSearch with data (visited'-'^ ?i;j) would also terminate with success as the last part of its' 
output. But because j < N , this cannot happen by Lemma 18-1 from which it follows that "W is disjoint 
from i?(visited(-') , ?i;j) . Thus W is also a probe walk for (visited(''+^), ?;) . 

By induction, "W is a probe walk for (visited'^^\ v) , so it must be that N ^ M . By the choice of A'' , there is 
then a probe walk W for (visited*^^\ ?i;jv) which ends in O . 
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Because w^f is part of the input to a daughter invocation of AugmentSearch, we have visited^^'(tt'_v) < iter; 
thus we can easily extend W (by one or two vertices, depending on whether wn is a neighbor of v or the 
predecessor in ^ of a neighbor of v) to form a probe walk W" for (visited^^^ z;) . Because W" will also be 
a probe walk for (visited, w) , it has length at most £: then W is strictly shorter than £ in length. By the 
induction hypothesis, the invocation of AugmentSearch on input data (visited*-^-*, wat) then halts, and returns 
the output value ® W , visited, success) , where W is a probe walk for (visited^^^ , Wn) ending in O , and 
where visited differs from visited^^^ only on a subset of R{v±s±ted^^\ wn) ■ We proceed by cases: 

• If t; is covered by a path of , v ^ I , and wat is the predecessor of u in then W = vwn'W' is a probe 
walk for (visited, v) ending in O . Note that A{J^ ® W) ~ A{J^ (BW ) \ {wjv — > v} ; then, the value which 
is returned as output on line 7 is (^ © W, visited, success) . 

• If Wn is not covered by a path of , then wn ^ v, and the walk W = vwnV^ is a probe walk for 
(visited, v) ending in O . Note that A{J^ © W) = A{f © W ) U {u — >■ wn} ; then, the value which is 
returned as output on line 16 is (^ © "W, visited, success) . 

• If neither of the previous two cases apply, it must be that wn is the predecessor in of some third vertex 

■u ~ w . Because wn is part of the input to a daughter invocation for (visited*^^\ , we know that 
visited(^^(M) < iter: then, W = wuwjvW is a probe walk for (visited,!;) ending in O . Note that 
A{J^ © W) = [A{J^ © W) \ {wn ^ u}] U — > zt} ; then, the value which is returned as output on line 23 
is (^©W, visited, success) . 

Finally, because i?(visited'^^\ wjv) ^ i?(visited, , and because visited'''^' differs from visited only on 
i?(visited(i), Wi) C i?(visited, u) for 1 < i < A/' , it follows that visited differs from visited only on a subset 
of i?(visited, w) , with visited(x) = iter on that subset. Thus, if the Lemma holds for pairs (visited, u) 
having probe walks of length less than £ ^ ending in O , it also holds for such pairs with probe walks ending 
in O of length £ + 1 . By induction, the Lemma then holds. 

To prove the Theorem, it then suffices to note that for a function visited : V{G) — > N with visited(a;) < iter 
for all X G V{G) , probe walks for (visited, i) are just proper alternating walks with respect to which start at 
i , in which case such a probe walk W which ends in O is a proper augmenting walk for . Then all the various 
parts of the Theorem follow from Lemmas 18-1 and 18-2 collectively. □ 

Run-time analysis. Because AugmentSearch marks each vertex v with visited(z;) <— iter when it visits v , 

each vertex is only visited once. At each vertex, each of the neighbors w ^ v are tested for if they fulfill the 
condition of line 3, or of lines 11, 12, and 18. Because computing prev , AddArc , and RemoveArc is constant-time 
for a collection of vertex-disjoint paths (or differing only slightly from one as described in the discussion 
on implementation details), the aniomit of work in an invocation to AugmentSearch for a vertex v £ V{G) is 
0{degv) , neglecting the work performed in descendant invocations. Summing over all vertices v G V{G) , the 
run-time of AugmentSearch is then 0(m) for an input as described in the statement of Theorem 18. 

3.1.3 An efficient algorithm for constructing a path cover for {G,I,0) 

Using AugmentSearch as a subroutine to build successively larger proper families of vertex-disjoint I - O paths. 
Algorithm 2 describes a straightforward subroutine which attempts to build a path cover for (G, /, O) . 

Corollary 19. Let {G,I,0) be a geometry with \I\ = \0\ : then BuildPathCover halts on input {G,I,0) . 
Furthermore, let a = BuildPathCover(G, 7, 0) . If a = fail, then {G,I,0) does not have a causal flow; 
otherwise, a is a path cover T for {G, I, O) . 

Proof — Suppose {G,I,0) has a causal flow: then it has a collection of fc |/| = \0\ vertex-disjoint / - 
O paths by Lemma 3. Then, by Corollary 15, for any proper collection of vertex-disjoint I - O paths with 
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Algorithm 2: BuildPathCover(G, I, O) — tries to build path cover for (G, J, O) 
Require: (G, /, O) is a geometry. 

1: let an empty collection of vertex-disjoint dipaths in G 

2: let visited : V{G) — > N be an array initially set to zero 

3: let iter ^ 

4: for alH e / do 
5: iter <— iter + 1 

6: {T, visited, status) ^ AugmentSearch(G, I, O, iter, visited, i) 
7: if status = fail then return fail 

8: end for 

9: if V{G) \ V{J^) = then 
10: return 

11; else 

12: return fail 

13: end if 



1^1 < k , there is a proper augmenting walk for !F starting at any i G I which is not covered by . For such a 
collection and vertex i , if visited(x) < iter for all x G V{G) , AugmentSearch(G, /, O, iter, visited, i) 
returns © W, visited, success) , where visited(x) ^ iter for all x € V{G) , and where W is a. proper 
augmenting walk for starting at i . Then, J^© T4^ is a proper collection of vertex-disjoint paths, covering i and 
the input vertices covered by , and with |JF0 W\ = |JF| -I- 1 . By induction, we may then show that at the 
end of the for loop starting at line 4, will be a family of vertex-disjoint I - O paths which covers all of J , in 
which case \!F\ = k . If all of the vertices of V{G) are covered by ^ , ^ is then a path cover for (G, /, O) , and 
BuildPathCover returns . Taking the contrapositive, if BuildPathCover(G, I, O) returns fail, then (G, 7, 0) 
has no path cover. 

Conversely, if BuildPathCover(G, /, O) returns fail, then either the condition of line 7 failed, or the condition 
of line 12 failed. If the former is true, then by Theorem 18 there were no proper augmenting walks for some 
proper collection of fewer than k disjoint I-O paths, in which case by Corollary 15 there is no such collection 
of size k , and thus no causal path cover for (G, /, O) . Otherwise, T is & maximum-size collection of disjoint 
paths from 7 to O , but is not a path cover for (G, /, O) ; then by Theorem 9, there again is no causal path cover 
for (G, 7, 0) . In either case, there is no causal flow for (G, 7, 0) by Theorem 8. The result then holds. □ 

Run-time analysis. BuildPathCover iterates through k = |7| input vertices as it increases the size of the 
collection of vertex-disjoint paths, invoking AugmentSearch for each one. The running time for this portion of the 
algorithm is then 0{km) . As this is larger than the time required to initialize visited or to determine if there 
is an element v € V(G) such that v ^ V{J^) , this dominates the asymptotic running time of BuildPathCover. 

3.2 Efficiently finding a causal order for a given successor function 

Given a path cover C for a geometry (G, I, O) , and in particular the successor function / of C , we are interested 
in determining if the natural pre-order =^ for / is a partial order, and constructing it if so. In this section, I 
present an efBcient algorithm to determine whether or not =^ is a partial order, by reduction to the transitive 
closure problem on digraphs. 

3.2.1 The Transitive Closure Problem 

Any binary relation R can be regarded as defining a digraph D with {x y) € ^(D) <^=> (xRy) . Chains 
of related elements can then be described by directed walks in the digraph D . This motivates the following 
definition: 
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Definition 20. Let / be a successor function for a geometry (G, /, O) : the influencing digraph 3 j is tlicn the 
directed graph witli vertices V"(3/) = V{G) , where {x ^ y) & A{3f) if one oiy = x , y = f{x) , or y ~ f(x) liold. 

The three types of arcs in Definition 20 correspond to the relations in Equation 2, whose transitive closure is the 
natural pre-order. Note that aside from self-loops x ^ x , the arcs in 3/ correspond directly to the two varieties 
of segments of influencing walks. (This is an alternative way of proving Lemma 6.) 




Figure 6: On the left: a geometry (G, /, O) with a path cover C . Arrows represent the action of the successor function / ; O'^ > /'^ 

of C . On the right: the corresponding influencing digraph Jf . Solid arrows represent arcs of the form x —* f{x) , and hollow arrows 
represent arcs a; — > y for y ~ f{x) . (Self-loops are omitted for clarity.) 

It is natural to also speak of transitive closures of binary relations in graph-theoretic terms, as follows: 

Definition 21. For a digraph D , the transitive closure of D is the digraph T with V{T) = V{D) , and such 
that (a; — > ?/) € A{T) if and only if there is a non-trivial^ directed walk from a; to in the digraph D . 

Thus, x ^ 2/ is an arc in the transitive closure of J/ and only ii x =4 y , ov equivalently iff there is an influencing 
walk for C from x to y in G. 

Transitive Closure Problem. Given a digraph D , determine it's transitive closure T . 

The Transitive Closure problem is known to be efliciently solvable. Algorithm 3 presents on solution, which is (a 
paraphrasing of) the pseudocode of Figure 3.8 from [13]. This algorithm is a simple modification of Tarjan's al- 
gorithm for finding strongly connected components of digraphs (equivalence classes of mutually reachable vertices 
using directed walks), which finds the transitive closure by determining the "descendants" of each x e V{D) : 

Desc(a;) = {y & V{D) | D contains a non-trivial directed walk from x to y} . (3) 

The following is an overview of Algorithm 3: interested readers may refer to [13] for a more complete analysis. 

• A di-connected component of D is an equivalence class of vertices which can be reached from each other by 
non-trivial directed walks in D . Tarjan's algorithm detects these components by performing a depth-flrst 
search which traverses arcs of D , and detecting when it has traversed a directed cycle in D . 

• A stack is used to keep track of vertices of the digraph have been visited, but whose di-connected compo- 
nent has not yet been completely determined. When the vertices belonging to a a given component are 
determined, we pop them off of the stack (line 14) and insert them into a set representing that component. 

• We say that v precedes w in the ordering of the stack if v is on the stack and w is not, or if v is lower on 
the stack than w is. Then, we may keep track of the "root" Root(w) of v , which is an upper bound on the 
stack-minimal vertex of the component containing v. At first, we set the root of v to itself, and we always 
ensure that Root(w) < v . 

Suppose we discover a descendant w of w such that Root(u>) ^ Root(u) ^ v . Then v is a descendant of 
Root(w;) , which is in a common component with w by definition. Because w is also a descondent of u , w 

''Note that if a vertex x has a loop a; — > a: (which are permitted in digraphs), then the directed walk a; ^ a: is a non-trivial walk. 
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Algorithm 3: Figure 3.8 of [13] — an algorithm for transitive closure of a digraph 

1: procedure SimpleTC(v) 

2: begin 

3: Root(w) ^ V ; Coinp(u) nil 

4: PUSH(i;, stack) 

5: Desc(u) ^ {w e V{D) \ (v ^ w) G A{D) } 

6: for all w such that (v w) G A{D) do 

7: if {w is not already visited) then SimpleTC(w) 

8: if Comp(w) = nil then Root(t;) ■*— min(Root(w), Root(w)) 

9: Desc(w) <— Desc(u) U Desc(t«) 

10; end for 

11: if Root(u) = V then 

12: create a new component C 

13: repeat 

14: let w <— POP(stack) 

15: Comp(w) ^ C 

16: insert w into the component C 

17: Desc(w) <— Desc(w) 

18: until w = V 

19: end if 

20: end 

21: procedure main 

22: begin 

23: let Stack ^ 

24: for all V e V{D) do 

25: if {v is not already visited) then SiinpleTC(t;) 

26: end for 

27: end 



must be in a common component with w . Then Root(ti)) is the smallest known vertex in that component: 
we update Root(u) <— Root(w) to improve the known minimum for v . 

• Because vertices are only allocated to a di-connected component after they are popped off the stack, we 
may test each of the descendants w of w to see if they have been allocated to a component, rather than 
testing if Root(w) ^ Root(ij) . If not, then w is in a common component with w , and we update Root(t;) 
to be the minimum of Root(f ) and Root(i(;) on line 8, as in the previous case. 

• If Root(t') = V on line 11, then v is the stack-minimal element of its' component: then any vertices higher 
than V on the stack will be in the same component as v . Conversely, because all descendants of v have 

been visited by that point, all of the vertices in the same component as v are still on the stack. Thus, we 
may pop them off the stack and allocate them to a component, until we have removed v off of the stack 
(lines 11 through 19). 

• As we determine the connected components of the digraph, we may maintain the sets of descendants of 
each vertex: if {v ^ w) G A{D) , then the descendants of w are all also descendants oi v , so we ensure 
that Desc(w) C Desc(v) (as on line 9). 

The above is performed for all vertices v G V{G) to obtain the transitive closure. 

Algorithm 3 is sufficient to build the natural pre-order =^ for a successor function / . However, the output does 
not indicate whether =^ is a partial order, and it performs work that is unnecessary if =^ is not actually a partial 
order. We may also take advantage of the availability of the path cover C which is given as input, which is not 
available in the more general Transitive Closure problem. Therefore, we are interested in adapting Algorithm 3 
to the application of finding a causal order. 
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3.2.2 Chain decompositions with respect to the path cover C 

Let C be a path cover for (G, /, O) with successor function / . The transitive closure of the influencing digraph 3f 
will often have high maximum degree: because the longest path in C has at least n/k vertices, and the end-point 
of this path will be at the terminus of arcs coming from every vertex on the path, the maximum in-degree of the 
transitive closure is at least n/k; and similarly for the maximum out-degree. In Algorithm 3, this implies that 
the set Desc(w) may become comparable to V{G) in size. In order to construct the arc-lists of the transitive 
closure reasonably efficiently, we want to reduce the effort required in determining the sets 0630(11) . 

A standard approach to this problem would be to find a chain decomposition [13] for 3/ , which is a collection of 
vertex-disjoint dipaths of J/ which cover all of J/ . By the definition of the influencing digraph, C itself is such 
decomposition of J/ . Then, using a chain decomposition with respect to C , we can efficiently represent Desc(a;) 
in terms of the first vertex y in each path of C such that y GT)esc{x) . 

Definition 22. Let C = {CPj}^.g^ be a parameterization of the paths of a path cover C for a geometry (G, /, O) , 
let / be the successor function of C , and let =^ be the natural pre-order for / . Then, for x G V{G) and j G K , 
the supremum supj{x) of x in is the minimum integer m e N, such that a; =^ y for all vertices y e VC^j) 
which are further than distance m from the initial vertex of . 

We may use the suprema of x in the paths of C to characterize the natural pre-order for / : 

Lemma 23. Let C = {^Pj } be a parameterization of the paths of a path cover C for a geometry (G, /, O) , 
let f be the successor function of C , let =4 be the natural pre-order for f , and let L : V{G) — > N map vertices 
X e V{G) to the distance of x from the initial point of the path of C which contains x . Then 

x^y <^=^ supj{x)^L{y) (4) 

holds for all x € V{G) and y € V{7'j) , for any j G K . 

Proof — Let x G V{G) , and fix G C. Let v G V^(J'j) be such that L{v) = sup^(a;) . By definition, if 
y G V{yj) and x =^ y , then L{y) ^ L{v) . Conversely, if y G VCPj) and L{y) = L{v) + h for h ^ Q , then 
y = f^{v) ; then x =4 v =4 y , and the result holds by transitivity. □ 

To determine the supremum function for all vertices, it will be helpful to be able to efficiently determine which 

path of C a given vertex belongs to and how far it is from the initial vertex for it's path. Algorithm 4 describes a 
simple procedure to do this, which also produces the successor function for the path cover C . (In the case where 
|/| = \0\, every path of C has an initial point in 7 ; we then take K = I to he the index set of the paths of C) 

3.2.3 Detecting vicious circuits with respect to C 

If the influencing digraph 3f contains non-trivial di-connected components, we know that there are closed in- 
fluencing walks — i.e. vicious circuits - for C in (G,I,0) . In that case. Theorem 8 together with Theorem 9 
imply that (G, /, O) has no causal flow, in which case we may as well abort. Recall that SimpleTC keeps track of 
di-connected components by allocating vertices to a component C after the elements of C have been completely 
determined. However, the state of being allocated into a component can be replaced in this analysis by a,ny 
status of the vertex which is changed after the descendants of a vertex have been determined; and this status 
may be used to determine if a vicious circuit has been found. 

Algorithm 5 is a simple procedure to initialize an array status over V{G) . A status of none will indicate that 

no descendants of the vertex have been determined (except itself), fixed will indicate that all descendants of 
the vertex have been determined, and pending will indicate that the descendants are in the course of being 
determined. Because output vertices have only themselves for descendants, their status is initialized to fixed; 
all other vertices are initialized with status(u) = none . At the same time. Algorithm 5 initializes a supremum 
function which represents only the relationships of each vertex to the ones following it on the same path. 
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Algorithm 4: GetChainDecomp(G, /, O, C) obtain the successor function / of C, and obtain functions 
describing tiie ciiain decomposition of the influencing digraph 3 f 

Require: {G,I,0) be a geometry with |7| = \0\ 
Require: C a path cover of {G, I, O) 

1: let P : V{G) — > / an array 

2: let L:V{G) — > N an array 

3: let / : O*^ — > I'^ an array 

4: for all i G / do 

5: let V ^ i , £ ^ 

6: while V do 

7: /(w) ^ next(C, w) 

8: P{v) ^ i ; L{v) ^ £ 

9: V ^ f{v) 

10: £^£+l 

11; end while 

12: P{v) ^ i ; L{v) ^ e 

13: end for 

14: return (/, P, L) 

Algorithm 5: InitStatus(G', /, O, P, L) — initialize the supremum function, and the status of each vertex 
Require: (G, 7, 0) is a geometry 

Require: P : V{G) — > / maps each x £ V{G) to i E I such that x is in the orbit of i under / 
Require: L : V{G) — > N maps each x G V{G) to ft, € N such that x = f^{P{x)) 

1: let sup : / X V{G) — > N an array 

2: let status : V{G) — > {none, pending, fixed} an array 

3: for all V £ V{G) do 

4: for all i G / do 

5: if i = P{v) then sup(i, v) <— L{v) 

6: else sup(i,i;) <— 

7: end for 

8: if ?; G O then status(t;) <— fixed 

9: else status(w) <— none 

10: end for 

11: return (sup, status) 



3.2.4 An efficient algorithm for computing the natural pre-order of / 

Algorithms 6 and 7 below represent a modified version of Algorithm 3, specialized to the application of computing 

the natural prc-ordcr for the successor function / of a path cover C . Rather than explicitly constructing the 
influencing digraph "if and traversing directed walks in Jj (as is done in Algorithm 3), we instead traverse 
influencing walks for C (characterized by its' successor function) in the graph G . 

Theorem 24. Let f be a successor function of a path cover C for a geometry (G, /, O) . Let P : V{G) — > / 
map vertices v to the initial point of the path of C that covers v , and let L : V{G) — > N map vertices v to the 
integer /i G N such that v — f^{P{v)) . Then ComputeSuprema halts on input {G, I,0, f, P, L) . Furthermore, 
let a = ComputeSuprema(G, /, O, /, P, L) . If cr = fail, then (G,I,0) does not have a causal flow; otherwise, 
(G, /, O) does have a causal flow, and a is a supremum function sup : I x V{G) — »■ N satisfying 

x4y <S=^ s\ip{P{y),x) ^ L{y) (5) 

for all x,y (z V{G) , where =^ is the natural pre-order for f . 

Proof — We will reduce the correctness of Algorithms 6 and 7 to that of Algorithm 3, where D = Jj is the 
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Algorithm 6: Traverseinf lWalk(G, /, O, f, sup, status, v) — compute the suprema of v and all of its' de- 
scendants, by traversing influencing walks from v 
Require: {G, I, O) is a geometry 

Require: / : O'^ — > I'^ is a successor function for [G, 7, 0) 
Require: sup : / x V{G) — > N 

Require: status : V{G) — > {none,pending, fixed} 
Require: u G O*^ 

1: status(w) ^ pending 

2: for all w = f{v) and for all w ~ f{v) do 

3: if w ^ V then 

4: if status(w) = none then (sup, status) <— Traverseinf lWalk(G, I, O, /, sup, status, w) 
5: if status(w) = pending then 
6: return (sup, status) 

7: else 

8: for alH e / do 

9: if sup(i,t;) > swp{i,iu) then sup(i,i)) <— swp{i,w) 

10: end for 

11: end if 
12: end if 
13: end for 

14: status(t!) <— fixed 
15: return (sup, status) 

Algorithm 7: ComputeSuprema(G, 7, O, /, P, L) — obtain the successor function / of C, and compute the 
natural pre-ordcr of / in the form of a supremum function and functions characterizing C 
Require: {G,I,0) is a geometry with |7| = \0\ and successor function / : — > 7'^ 
Require: P : V{G) — > 7 maps each x € V{G) to i e I such that x is in the orbit of i under / 
Require: L : V{G) — > N maps each x e V(G) to /i e N such that x = f^{P{x)) 

1: let (sup, status) <— InitStatus(G, 7, O, P, i) 

2: for all V do 

3: if status(w) = none then (sup, status) <— Traverseinf lWalk(G, 7, 0, /, sup, status, v) 

4: if status(w) = pending then return fail 

5: end for 

6: return sup 



digraph provided as the the input of the main procedure. Throughout, =^ denotes the natural pre-order of / . 

For distinct vertices w, w G y{G) , because (w ^ w) G A{3 f) if and only if either w = f{v) oi w ^ f{v) , wc may 
replace the iterator limits "w such that {v ^ w) € A{D)" of the for loop starting on line 6 of Algorithm 3 with 
a loop iterating over w = f{v) and w ^ f{v) : this is what we have on line 2 of Traverseinf IWalk. 

At line 8 of Algorithm 3, if Comp(u') = nil, we infer that v and w arc in a common di-connected component of 
the digraph 3 j : this implies that v =<; 7u and w ^ v . If v ^ w . this implies that =<; is not antisymmetric, and 
thus not a partial order; by Theorem 8, C is then not a causal path cover. We proceed by cases: 

• If =^ is antisymmetric, then the influencing digraph is acyclic, in which case 3 y has only trivial di-connected 
components. In this case, the following changes preserve the functionality of Algorithm 3: 

— In the case that w = t; in the for loop, all the operations performed are superfluous, in which case we 
may embed lines 7 through 9 in an if statement conditioned o'n.w i^v . 

— Because each vertex is the only vertex in its' component when 3/ is acyclic, we may replace lines 11 

through 19 of SimpleTC with a line setting Comp(u) to an arbitrary non-nil value, which in this case 
may be interpreted as allocating the vertex v to its' di-connected component (i.e. the singleton {y\). 
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Also, the condition of line 8 is never satisfied in a call to SimpleTC(w) . Then, we may replace the 
conditional code with an arbitrary statement, e.g. a command to abort the procedure. 

— After the above replacement, the value of stack is not used within the procedure call SimpleTC(v) , 
and has the same value after the procedure call to SimpleTC(f ) as it docs before the call. Then, 
stack is superfluous to the performance of the algorithm. Similarly, the value of Root(t;) is not 
affected except to initialize it. We may then eliminate all references to either one. 

— The value of Coinp(w) is only tested to determine whether or not it is nil , so we may replace the array 
Comp with status . and its' possible states of being nil or non-nil with the states of being pending 
and non-pending. We define the two values none and fixed to represent being non-pending and also 
having not yet been visited, and being non-pending and having been visited, respectively. 

— Using the array sup to implicitly represent the sets of descendants, we may replace the union performed 
on line 9 with code which sets sup(i, v) to the minimum of sup(i, v) and sup(i, w) for each i & I . Note 
also that, because x is a descendant of u in 3/ iff u = a; or a; is a descendant of w , we may remove 
the initialization of Desc(w) on line 5 of Algorithm 3 if we initialize sup for each vertex so that it 
represents each vertex as a descendant of itself (for instance, in the main procedure, which is replaced 
by ComputeSuprema). 

By performing the substitutions described above, we can easily sec that Traverseinf IWalk together with 
ComputeSuprema is equivalent to Algorithm 3 when =<! is anti-symmetric. Then, a ^ fail because line 4 
of ComputeSuprema is never evaluated; we then have a = sup as in Equation 5, from the correctness of 
Algorithm 3. 

• If =^ is not antisymmetric, then there are distinct vertices x,y € V{G) such that x ^ y ^ x , in which case 
X and y are in a non-trivial component in J/ . Then, the for loop of BuildCausalOrder will eventually 

encounter a vertex v of which x and y are descendants. 

In the depth-first traversal of influencing walks performed in Traverseinf lWalk(G, /, O, f, sup, status, v) , 
eventually a directed cycle containing both x and y will be discovered. Without loss of generality, assume 
that the depth-first traversal starting from v visits x before y : then, the depth-first traversal will eventually 
uncover a walk of the form 

Then in the procedure call Traverseinf lWalk(G', /, O, /, sup, status, y') , line 5 will find status(x) = 
pending, as line 1 of the procedure call Traverseinf lWalk(G, J, O, /, sup, status, a;) has been executed 
while line 14 has not. Then, the procedure aborts by returning TZ without first changing the status of 
status(j/') from pending. 

It is clear that if w' depends on w , and if Traverseinf lWalk(G', /, O, /, sup, status, u/) aborts with 
status(u'') = pending during a procedure call Traverseinf lWalk(G, /, O, /, sup, status, w) , then the 
latter will also abort with status(t«) = pending. By induction, we may then show that for v G V{G) for 
which X and y arc descendants, Traverseinf lWalk(G, /, O, f, sup, status, v) will abort with status(u) = 
pending in the for loop in ComputeSuprema. 

By the analysis of the case where ^ is antisymmetric, the status status(u) = pending will only occur at 
line 4 of ComputeSuprema if =^ is not antisymmetric. If this occurs, a = fail ; as well, no causal path cover 
exists for (G, /, O) by Theorem 9, and thus no flow exists for (G, I, O) by Theorem 8. 

Thus, a ^ fail iff =^ is a partial order; and when this occurs, by reduction to Algorithm 3, sup corresponds to 
the natural pre-order =<; in the sense of Equation 5. □ 

Run-time analysis. We may analyze the run-time of Algorithm 7 as follows. Let n = \V{G)\ , m = \E{G)\ , 
k= \I\ = |0| , and d be the maximum degree of G . The time required to execute the for all loop starting on line 8 
of Traverseinf IWalk is 0{k) ; then, aside from the work done in recursive invocations to Traverseinf IWalk, 
the time required to perform an invocation of Traverseinf IWalk for a vertex v is 0{kdegf{v)) . Because the 
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first invocation of Traverseinf IWalk for a vertex v will change status(t;) to something other than none, which 
prevents any further invocations for v , Traverseinf IWalk will only be called once for any given vertex in the 
course of Algorithm 7. Then, summing over all vertices v S V{G) , the amount of time required to perform the 
for all loop starting on line 2 of ComputeSuprema is 0{km) . The time required by InitStatus to initialize sup 
and status is 0{kn) ; then, the overall running time of Algorithm 7 is 0{km) . 

3.2.5 A slightly more efficient algorithm for finding a causal order for / 

If C is a causal path cover, it also is possible to find a causal order compatible with / which differs from the 
natural pre-order for /, or determine that none exists, by recursively assigning integer "level" values to vertices 
rather than building the set of descendants. For example, one may construct a function A : V{G) — > N satisfying 

A(a;) = , if there are no influencing walks for C ending at x; 

\{x) = l+max{A(y) \ x = f{y) or a; ~ /(?/)}, otherwise. 

Note that the set S{x) of vertices y such that x = f{y) or x ~ f{y) are the initial points for any influencing walk 
for C with one segment which ends at x. By constructing the predecessor function g — f^^ of C rather than 
the functions P and L in Algorithm 4, we can easily find all elements of S{x) in G by visiting g(z) for z — x or 
2; ~ a; . Then, such a level function can be constructed by a Tarjan style algorithm similar to Algorithm 6, using 
the status array in the same way, but traversing the arcs of the influencing digraph J/ in the opposite direction 
as Traverseinf IWalk. We may then define a; y [a; = y] V [X{x) < X{y)] . 

It is easy to see that the resulting partial order resulting would have the same maximum-chain length as the 
natural pre-order =^ : any maximal chain in =^ is a list of the end-points of consecutive segments in an influencing 
walk for C , which will be a maximum chain in \^\. However, also contains relationships between vertices with 
no clear relation in the influencing digraph 3 / , because it suffices for two vertices to be on different "levels" for 
them to be comparable. 

Such a causal order can actually be constructed in 0{m) time, because the algorithm to construct it consists 

essentially of just a depth-first traversal with operations taking only constant time being done at each step. We 
have instead presented the above algorithm because the extra time required to obtain the coarsest compatible 
causal order for / will not affect the asymptotic run time of the complete algorithm for finding a flow, because 
of the immediate reduction to the well-studied problem of transitive closure, and in the interest of describing an 
algorithm to construct the natural pre-order for / (being the coarsest compatible causal order for /). 

3.3 The complete algorithm 

We now describe the complete algorithm to produce a flow for a geometry (G, /, O) , using Algorithms 2 and 7. 



Algorithm 8: FindFlow(G, I, O) — try to find a flow for (G, I, O) 
Require: {G,I,0) is a geometry with |7| = \0\ 

1: let T ^ BuildPathFamily(G,/, O) 
2: if T = fail then return fail 

3: let (/, P, L) <- GetChainDecomp(G, /, O, r) 

4: let a ^ ComputeSupreina(G, I, O, /, P, L) 

5: if (7 = fail then 
6: return fail 

7: else 

8: return (/, P, L, a) 

9: end if 
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Corollary 25. Let {G,I,0) be a geom,etry with \I\ — \0\ . Then FindFlow halts on input {G,I,0) . Further- 
more, z/ FindFlow(G, /, O) = fail, then {G,I,0) does not have a causal flow; otherwise, FindFlow(G, /, O) = 
(/,P,Z/, sup) , and {f,=4) is a causal flow, where =4 is characterized by 

x^y ^ snp{P{y),x) ^ L{y); (7) 

and is the natural pre- order for f . 

Proof — By Corollary 19, a causal path cover exists for (G, /, O) only if BuildPathFainily(G, /, O) sets 
T ^ fail on line 1; thus if r = fail , (G, I, O) has no causal flow by Theorem 8. Otherwise, r is a path cover. If 
BuildCausalOrder sets a = fail on line 4, (G, /, O) has no causal flow by Theorem 24. Otherwise, the relation 
=^ characterized by Equation 7 is the natural pre-order for / and a causal order, in which case (/, is a causal 
flow. □ 

Run-time analysis. Because r 7^ fail at line 1 implies that r is a path cover, GetChainDecomp visits 
each vertex v e V{G) once to assign values for P{v) , L{v) , and possibly f{v) in the case that v G O"^ . 
Then, its' running time is 0{n) . The running time of FindFlow is then dominated by BuildPathCover and 
ComputeSuprema, each of which take time 0{km) . 

4 Potential Improvements 

This paper has described eSicient algorithms for finding flows, with the aim of not requiring prior knowledge 
of graph-theoretic algorithms in the presentation. This constraint has led to choices in how to present the 
algorithms which may make them less efficient (in practical terms) than may be achievable by the state of the 
art; and no significant analysis of the graphs themselves have been performed. Here, I discuss issues which may 
allow an improvement on the analysis of this article. 

4.1 Better algorithms for finding path covers 

For network-fiow problems (the usual tools used for solving questions of maximum-size collections of paths in 
graphs), there is a rich body of experimental results for eflicient algorithms. However, there seems to be very 
little discussion in the literature of the special case where all edge capacities are equal to 1 , which is relevant 
to the problem of finding maximum collections of vertex-disjoint I - O paths. It is difficult to determine, in 
this case, whether there is a significant difference in the performance of various algorithms. Although it is 
less efficient than other algorithms for general network fiow problems, the most obvious choice of network flow 
algorithm for finding a maximum family of vertex-disjoint I O paths is the Ford-Fulkerson algorithm, which 
has an asymptotic running time 0{km) . This running time is identical to Algorithm 2: this should not be 
surprising, as Algorithm 1 essentially implements a depth-first variation of the Ford-Fulkerson algorithm for 
finding an augmenting flow. 

A more thorough investigation of network flows may yield an improved algorithm for flnding a path cover for 
(G, /, O) , which (when coupled with the faster algorithm for finding a minimum-depth causal order) would yield 
a faster algorithm for flnding causal flows. 

4.2 Extremal results 

Consider all the ways we can add edges between n vertices to get a geometry with k output vertices and a causal 
flow. Just to achieve a path cover, we require n — k edges; this lower bound is tight, as graph consisting of just 
k vertex-disjoint paths on n vertices has this many edges, and the paths represent a causal path cover of that 
graph. The more interesting question is of how many edges are required to force a graph to not have any causal 
path covers. 



23 



Let n be the residue of n modulo k . Consider a collection of k paths {Pj}^.gjj,j , given by Pj = p*"' ■ • .p*/"''''! ^' 
for j < n, and Pj = p^p ■ ■ -pj-"^^^'^^ for j ^ n. Then, let G be the graph defined by adding the edges p^^^^p^ 
for all a and h ^ j where these vertices are well-defined, and p^^'p/T^^' for all and h < j where these vertices 
are well-defined. We may identify the initial point of the paths Pj as elements of I and end-points as elements 
of O : then, let M{n, k) denote the geometry (G, /, O) constructed in this way. 

The geometry M{n,k) has the obvious successor function given by /{Pj"^) = Pj'^^^ for all j and a where both 
vertices are defined. Then, consider the natural pre-order for / : 

(«). we obviously have =^ Pj'' for a ^ 6 , for every j e [k] ; 

(n). from the edges p^"^ p^"' , we obtain p^""^' =^ p^"' for all h, j e [k] and all a > ; and 

(iii). from the edges P/j*"*"'^ pf^ for h < j , we obtain p^''^ =^ p^''' and pf~^^ =4 Ph^^'' ■ (Note that the second of 
these two constraints is redundant, as Pj"'^^ =4 P^^ =4 P^^^' is implied by the above two cases.) 

Then, the natural pre-order ^ on M(n, k) is closely related to the lexicographical order on ordered pairs: pj^"' and 
Pj"^ are incomparable if they are both endpoints of their respective paths P/j and Pj , and otherwise p^"* =^ Pj"^ if 
and only if either a < b, or a = b and h ^ j . This is clearly a partial order, so M{n, k) has a causal flow: and 
it has kn — {^%^) edges in total. 

I conjecture that this is the maximum number of edges that a geometry on n vertices with k output vertices can 
have. If this can be proven, we can determine that certain geometries have no flows just by counting their edges; 
the upper bounds of this paper can then be improved to 0{k^n) . 

5 Open Problems 

To conclude, I re-iterate the open problems presented in [7] . 

1. The general case. When |/| > |0| , it is easy to see that a causal flow cannot exist, because no successor 
function / may be deflned. This leaves the case where |/| < |0| • If ^ = \0\ — \I\ , we may test sets 
dl C 7*^ with \dl\ = 5 to see if the geometry (G, / U dl,0) has a causal flow: doing this yields an 
0{kmn^) algorithm for finding a causal fiow for (G, /, O) . Is there an algorithm for finding causal flows in 
an arbitrary geometry with |7| ^ \0\ , whose run- time is also polynomial in (5 = |0| — |7| ? 

2. Graphs without designated inputs/outputs. Quantum computations in the one-way model may be 

performed by composing three patterns: one pattern to prepare an appropriate quantum state, a pattern 
to apply a unitary that state (in the vein that we have been considering in this article), and a flnal pattern 
which measures the resulting state in an appropriate basis. The composite pattern has no input or output 
qubits, and so has only the measurement signals as an output. The result of the computation would then 
be determined from the parity of a subset of the measurement signals. 

Given a graph without any designated input or output vertices, what constraints are necessary to allow a 
structure similar to a causal flow to be found, which would guarantee that deterministic n qubit operations 
in the sense of [5] can be performed in the one-way measurement model with the entanglement graph G ? 

3. Ruling out the presence of causal flows with only partial information about G . Arc there 
graphs G where it is possible to rule out the presence of a flow for (G, /, O) from a proper sub-graph of 
G , or given n = \V{G)\ , m = \E{G)\ , and k = \I\ = \0\ ? (This question obviously includes the extremal 
problem asked earlier.) 

4. Relaxing the causal flow conditions for Pauli measurements. Suppose that, in addition to 7 and 
O , we know which qubits are to be measured in the X axis and which are to be measured in the Y axis 
(corresponding to measurement angles and 7r/2 respectively). These qubits can always be measured first 

in a pattern, by absorbing byproduct operations on those qubits and performing signal shifting. However, 
the analysis of patterns in terms of causal fiows does not take this into account, as it is independent of 
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measurement angles. Is it possible to develop a natural analogue for causal flows which represents these 
qubits as minimal in the corresponding causal order, which may be efficiently found for geometries with 
|/| = \0\ or |/| < \0\ generally? 



The results of this article were inspired by the similarity between of the characterization in terms of causal flows, 
with aspects of graph theory related to Menger's Theorem in general, and the relationship between influencing 
walks and alternating walks in particular. Investigation into open questions involving efficient construction of 
causal flows or relaxations of them may benefit from additional investigation of this link. 
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